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WHY DO THE BRAKING INDICES OF PULSARS SPAN A RANGE OF MORE THAN 100 MILLIONS? 
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Draft version November 27, 2012 

ABSTRACT 

Here we report that the observed braking indices of the 366 pulsars in the sample of Hobbs et al. 
range from about —10 s to about +10 8 and are significantly correlated with their characteristic ages. 
Using the model of magnetic field evolution we developed previously based on the same data, we 
derived an analytical expression for the braking index, which agrees with all the observed statistical 
properties of the braking indices of the pulsars in the sample of Hobbs et al. Our model is, however, 
incompatible with the previous interpretation that magnetic field growth is responsible for the small 
values of braking indices (< 3) observed for "baby" pulsars with characteristic ages of less than 
2 x 10 3 yr. We find that the "instantaneous" braking index of a pulsar may be different from the 
"averaged" braking index obtained from fitting the data over a certain time span. The close match 
between our model-predicted "instantaneous" braking indices and the observed "averaged" braking 
indices suggests that the time spans used previously are usually smaller than or comparable to their 
magnetic field oscillation periods. Our model can be tested with the existing data, by calculating 
the braking index as a function of the time span for each pulsar. In doing so, one can obtain for 
each pulsar all the parameters in our magnetic field evolution model, and may be able to improve the 
sensitivity of using pulsars to detect gravitational waves. 

Subject headings: magnetic fields - methods: statistical - pulsars: general - stars: neutron 



1. INTRODUCTION 
Assuming the braking law of a pulsar as 



(1) 



where v is its spin frequency and nb is called its braking 
index, Manchester & Taylor (1977) found that 



suggests that < 3 indicates magnetic field growth of a 
pulsar, and vice versa, since v < and K > 0. This can 
be seen more clearly from (Zhang & Xie 2012b; hereafter 
Paper I), 

(5) 



v 6 nb v 



nb = vvjii 2 , 



(2) 



if K = 0. For the standard magnetic dipole radiation 
model with constant magnetic field (K = 0), rib = 3. 
Therefore nb ^ 3 indicates some deviation from the 
standard magnetic dipole radiation model with constant 
magnetic field. 

Blandford & Romani (1988) re-formulated the braking 
law of a pulsar as, 



(3) 



This means that the standard magnetic dipole radiation 
is responsible for the instantaneous spin-down of a pul- 
sar, but the braking torque determined by K{t) may be 
variable. In this formulation, nb ^ 3 does not indicate 
deviation from the standard magnetic dipole radiation 
model, but means only that K(t) is time dependent. 
From Equation ^ one can obtain, 



Kv 
A v 



(4) 



Assuming that magnetic field evolution is responsible for 
the variation of K(t), we have K = AB(t) 2 , in which A 
is a constant and B{t) is the time variable dipole mag- 
netic field strength of a pulsar. The above equation then 
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Many authors have thus used the observed nb < 3 of 
some very young pulsars to infer their possible magnetic 
field increase (e.g., Blandford & Romani 1988; Chan- 
mugam & Sang 1989; Manchester & Peterson 1989; Lyne 
et al. 1993, 1996; Johnston & Galloway 1999; Wu et al. 
2003; Lyne 2004; Chen & Li 2006; Livingstone et al. 
2006, 2007; Middleditch et al. 2006; Yue et al. 2007; 
Espinoza et al. 2011). Magnetospheric activities or fall- 
back accretion may also be able to produce nb < 3 (e.g. 
Menou et al. 2001; Alpar et al. 2001; Xu & Qiao 2001; 
Wu et al. 2003; Chen & Li 2006; Yue et al. 2007). This 
can also be explained by another model based on a de- 
crease in effective moment of inertia due to an increase 
in the fraction of the stellar core that becomes superfluid 
as the star cools via neutrino emission (Ho & Andersson 
2012). 

However, the above studies have not investigated the 
statistical properties of the braking indices of pulsars. 
Recently, Magalhaes et al. (2012) studied the observed 
braking indices of several very young pulsars statistically. 
They built a model of a pulsar's braking law to explain 
the ranges of these observed braking indices, and then 
predicted the possible values of the braking indices of 
several other very young pulsars withusing only their 
measured v and v. Testing these predictions and ap- 
plying the model to more pulsars for observational tests 
may shed light on our further understanding of pulsars. 

In Paper I, we have shown that the observed nb, for the 
large sample of pulsars reported by Hobbs et al. (2010, 
hereafter H2010), not only deviates from nb = 3, but also 
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Fig. 1. — (a) Distributions of the observed braking indices 
("obs = vv/i> 2 ). (b) Ratios of the observed and predicted braking 
indices (n pre ) from Equation (0 (equation (11) in Magalhaes et 
al. 2012). 



ranges from around — 10 8 to around +10 8 , i.e., spanning 
a range of more than 100 millions, which has not been 
explained so far. In this work, we will first apply the 
model of Magalhaes et al. (2012) to this sample, and 
then search for possible correlations of the observed rib 
of pulsars in the sample of H2010. We find a significant 
correlation between rib and the characteristic ages of pul- 
sars, which is explained satisfactorily with the magnetic 
field evolution model we developed in Paper I. Finally we 
will make a testable prediction of our model and discuss 
its physical implications. 

2. APPLICABILITY OF THE MODEL OF MAGALHAES ET 
AL. (2012) TO THE PULSARS IN H2010 

Magalhaes et al. (2012) used equation ((T|) to predict, 



iog(\n\sye 

bg(fi) 



(6) 



where D, = 2ttv, S = 2.3 x 10 20 Hz 1 / 2 G and is assumed to 
be the same for all pulsars, and 1.5 x 10 13 < £ < 34 x 10 13 
in units of Hz( 3 ~" b )/ 2 G determines the range of n prc 
for each pulsar. Magalhaes et al. (2012) showed that 
Equation © describes the ranges of rib observed for sev- 
eral very young pulsars very well. A natural question 
arises: can this model also explain the observed rib for 
other older pulsars of the sample of H2010? In Figured] 
we plot the distributions of the observed braking indices 
(^obs = vv/v 2 ) of the pulsars of the sample of H2010 
and the ratios between n b s and n prc . n bs spans a range 
from around — 10 8 to around +10 8 , as shown in Paper I. 
However, n a bs/n pTe spans a similar range and has almost 
the same numbers of positive and negative values. This 
means that the model of Magalhaes et al. (2012) cannot 
explain both the magnitudes and signs of the braking in- 
dices of the pulsars of the sample of H2010. Therefore an 
alternative model is needed to account for the observed 
braking indices of these pulsars. 

3. OBSERVED CORRELATIONS OF BRAKING INDICES 

In order to develop a model to account for the observed 
braking indices of the pulsars of the sample of H2010, we 
first explore possible correlations of the observed braking 
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Fig. 2. — Observed correlations of the braking index rib with 
(from top to bottom) u, u, u and t c . 



index rib — vvjv 1 with various spin-down parameters 
and their combinations, as shown in Figure [5J No or 
only weak correlation is found either between rib and v 
or between rib and v. Strong correlations are found both 
between rib and i> and between rib and the characteristic 
age r c = — the latter correlation is very significant 
and is the focus of study in the rest of this work. 

The rib ~ t c correlation shown in the bottom panel 
of Figure [5] reveals three phenomena. (1) For all young 
pulsars with 10 4 < r c < 5 x 10 4 yr, 10 < rib < 40 and 
rib is not correlated with r c . (2) Excluding these pulsars, 
the overall correlation is almost linear over about eight 
orders of magnitude for both rib and t c . (3) The overall 
correlation is almost the same for positive or negative rib- 

4. A NEW ANALYTICAL MODEL OF BRAKING INDEX 

In Paper I, we have shown that the observed corre- 
lations between v and other observables or parameters 
can be well explained by assuming that the braking law 
takes the form of Equation ([3]) and the evolution of the 
dipole magnetic field of a pulsar consists of a long-term 
decay modulated by one or multiple oscillation compo- 
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nents, i.e., 

B(t) = Bo(-t-T a (l + ksm(0 + 27r^)), (7) 

where to and -Bo are the initial decay time and dipolc 
magnetic field strength, a is the decay index, and k, cp 
and T are the oscillation amplitude, initial oscillation 
phase, and oscillation period, respectively. In Paper I we 
have shown that a > 0.5, but its upper limit cannot be 
well constrained with the existing data. In this work, we 
simply take a = 0.5 for simplicity; the conclusions of this 
work will not change for larger values of a. 

From Equation (J3]) , we can solve for the spin evolution 
of a pulsar, 

v- 2 = v^ 2 + 2 [ K(t)dt. (8) 

J to 

With K = AB(t) 2 and Equation 0, we have 

t c = -— = 1 n sa tin — , (9) 
2v 2Kv 2 t W 

where a = 0.5 and are assumed. It is noticeable 

that the oscillation term in Equation is not impor- 
tant in the above relation, since k <C 1. However the 
decay component of Equation J7]) makes the real age of 
a pulsar significantly shorter than its characteristics age, 
since normally t/to 1 (Zhang & Xie, 2012a). 

From Equation (|SJ we can easily obtain v(t) and v(t). 
Inserting v(t), v(t) and i)(t) into Equation ([2J, we have, 

n h = 3 + lni-(2-4/tC(t)) = 3+^(2-4/tC(*)), (10) 
to t 

where / = ^ and C(t) = cos(0 + 2tt-^). As we have 
shown in Paper I, the oscillatory term can be ignored 
in determining i), so young pulsars with r c < 10 5 have 
v > 0. Ignoring the oscillatory term in Equation ([T0|). 
we have 

n b » 3 + 2— = 3 + 21n— , (11) 
t t 

which for young pulsars gives rib > (since t > to) and 
rib is not well correlated with r c , in agreement with the 
data shown in the lower left corner of the bottom panel in 
Figured! For old pulsars, the oscillatory term dominates 
and thus Equation (fTUf gives, 

n b « ±4r c /, (12) 

since C(t) most likely takes values around +1 or —1. This 
immediately suggests that nb should be linearly corre- 
lated with r c and have equal probabilities to be either 
positive or negative, again in agreement with the corre- 
lation shown in the bottom panel of Figure [21 

In Figure [31 we plot the analytically calculated n b as a 
function of r c using Equations (J9)) and (TlOl) for k = 10~ 4 
and 10 -5 , respectively; we also assume to = 0.5 yr and 
T = 10 yr, the latter is the typical time scale of the 
quasi-periodic features in the power spectra of the tim- 
ing residuals reported in H2010. The upper boundary in 
each panel is obtained when C(t) = ±1 and the other 
values of are obtained when — 1 < C(t) < 1. Clearly 
the observed correlation shown in Figure [2] is well repro- 
duced. In particular, the lack of nb < for young pulsars 
with r c is predicted with k < 10~ 4 , which in our model 
is caused by their long term magnetic field decay. 



Equation (|10[) can be rewritten as, 

n h /r c = 3/r c + 2/t - 4/C(i) « 3/r c + 2/t ± 4/. (13) 

In Figure [4l we show the observed correlation between 
Ub/Tc and r c , overplotted with the analytical results of 
Equation {BJ with T = 10 yr and k = 10~ 3 , 10" 4 and 
10~ 5 , respectively. Once again, the analytical results 
agree quite well with the data quite well. 

5. ON "BABY" PULSARS AND YOUNG PULSARS 

As discussed above, < 3 has been well established 
for all very young pulsars with r c < 2 x 10 3 yr (Living- 
stone et al. 2007), which we call "baby" pulsars here for 
convenience. In Figure [5] we show a comparison between 
all the measured braking indices of "baby" pulsars and 
those of young pulsars with 10 < r c < 10 5 yr. From 
Equation (0, we get K > or K < for "baby" pul- 
sars or young pulsars, respectively. If K ^ is due to 
B 7^ 0, we are then led to the scenario that the magnetic 
fields of "baby" pulsars grow, but those of young pul- 
sars decay. If these "baby" pulsars with r c < 2 x 10 3 yr 
evolve to become young pulsars, this must happen when 
the time of transition from the growth phase to the de- 
cay phase is longer than the ages of these "baby" pul- 
sars, i.e., ^Q aby > 10 3 yr. For a typical young pulsar 
with r c ~2x 10 4 yr and ~ 20, Equation ([Pi]) gives 
t ~ 2.4 x 10 3 and ^ oung ~ 0.5 yr. The mismatch be- 
tween tp aby and £g 0uns suggests that significantly differ- 
ent mechanisms are responsible for the observed braking 
indices for "baby" and young pulsars, respectively. In 
other words, our model that nb ^ 3 is due to B ^ 
can only be applied to young (and old pulsars), but not 
to "baby" pulsars. Indeed, we may expect significant 
fall-back accretion or neutrino cooling only for "baby" 
pulsars; these mechanisms have previously also been sug- 
gested to be responsible for the observed nb < 3 of these 
"baby" pulsars. In either case, the additional torque 
of accretion or the reduced moment of inertia due to 
an increased fraction of superfluid stellar core can re- 
sult in nb < 3. Stronger magnetospheric activities are 
also expected from "baby" pulsars, whose magnetic fields 
should be stronger than those of young pulsars. 

6. MODEL PREDICTION: THE EFFECT OF TIME SPAN 

All the above equations for nb give nb as a function of 
t, i.e., the calculated nb is in fact a function of time for 
a given pulsar, as shown in the left panel of Figure [6j in 
which the horizontal axis "Time" is the calendar time. 
We call nb calculated this way the "instantaneous" brak- 
ing index. However, in analyzing the observed timing 
data of a pulsar, one usually fits the data on time of 
arrival (TOA) over a certain time span to a Taylor ex- 
pansion to third order: 

*(t) = $o + Mt-to) + li<o(t-to) 2 + lMt-to) 3 , (14) 
2 o 

where <I>(t) is the phase of TOA of the observed pulses, 
and $o , vo, an d vo are the values of these parameters 
at to, to be determined from the fitting, nb calculated 
from vo, vo and vq is thus not exactly the same as the 
"instantaneous" braking index. We call nb calculated 
this way over a certain time span the "averaged" braking 
index. 
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Fig. 3. — Analytically calculated n.v, as a function of r c using Equation lllOU , for different values of k. 
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Fig. 4. — Correlation between n^,/r c and t c . The different curves show the analytically calculated model predictions using Equation (1 1 3 I t 
for different values of k. 
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Fig. 5. — Comparison between the measured braking indices of 
"baby" pulsars (in the lower left corner; data from Livingstone et 
al. 2007) and those of young pulsars (in the upper right corner; 
data from H2010). 



In the right panel of Figure [BJ we show the simulated 
result for the "averaged" braking index as a function of 



time span, In this simulation, we first produce a series of 
TOAs using Equations ([U and (JTJ), and then use Equa- 
tion (I14[) to obtain Uq, vq and v§ for different lengths 
of time span. It can be seen that the "averaged" 
is close to the "instantaneous" one when the time span 
is shorter than T, which is the oscillation period of the 
magnetic field. The close match between our model pre- 
dicted "instantaneous" and the observed "averaged" 
rib , as shown in Figure [4] suggests that the time spans 
used in the H2010 sample are usually smaller than T. 

For some pulsars the observation history may be longer 
than T and one can thus test the prediction of Figure [6] 
with the existing data. In doing so, we can also obtain 
both / and T for a pulsar, thus allowing a direct test of 
our model for a single pulsar. We can in principle then 
include the model of magnetic field evolution for each 
pulsar in modeling its long term timing data, in order 
to remove the red noise in its timing residuals, which 
may potentially be the limiting factor to the sensitivity 
in detecting gravitational waves with pulsars. 

7. SUMMARY 
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Fig. 6. — Left: "instantaneous" braking index as a function 
of time. Right: "averaged" braking index as a function of the 
time span of the fitting. T = 15 yr is used in both cases. 



5. We find that the "instantaneous" braking index of 
a pulsar may be different from the "averaged" brak- 
ing index obtained from data. The close match be- 
tween our model-predicted "instantaneous" brak- 
ing indices and the observed "averaged" braking 
indices, as shown in Figure|4]suggests that the time 
spans used in the H2010 sample are usually smaller 
than or comparable to their magnetic field oscilla- 
tion periods. 



The results and conclusions of this work can be sum- 
marized as follows. 

1. The observed braking indices of the pulsars in the 
H2010 sample, which span a range of more than 
100 millions, are completely different from the pre- 
dictions of the model of Magalhaes et al. (2012), 
as shown in Figure [T] 

2. A significant correlation between the braking in- 
dices and characteristic ages of pulsars in the H2010 
sample is found over about eight orders of magni- 
tude, as shown in the bottom panel of Figure [21 

3. Based on the magnetic field evolution model we de- 
veloped previously (Paper I), an analytical expres- 
sion for the braking index is derived, as given by 
Equation (fT0|) . The analytically calculated correla- 
tion between the braking index and characteristic 
age re-produces the observed correlation well, as 
shown in Figured] 

4. Our model is incompatible with the previous in- 
terpretation that the magnetic field growth is re- 
sponsible for the observed small values of braking 
indices (< 3) of "baby" pulsars with characteristic 
ages of less than 2 x 10 3 yr. 



6. Our model can be tested with the existing data, 
by calculating the braking index as a function of 
the time span for each pulsar. In doing so, one 
can obtain all the parameters for each pulsar in our 
magnetic field evolution model. This is particularly 
important if one wants to use the long term timing 
data of pulsars to detect gravitational waves. 
The anonymous referee is thanked for a valuable sug- 
gestion which helped to clarify one important point 
in the revised manuscript. S.-N.Z. acknowledges par- 
tial funding support by 973 Program of China under 
grant 2009CB824800, and by the National Natural Sci- 
ence Foundation of China under grant Nos. 11133002, 
10821061, and 10725313. 

Note added. After our Paper I was published and this 
manuscript was submitted, we noticed that Pons et al. 
(2012) also proposed a phenomenologically similar model 
of magnetic field oscillations, which was also used to ex- 
plain the observed braking indices of older pulsars. In 
their model the oscillations are identified as due to the 
Hall drift effect and have time scales of (10 6 -10 8 )/f? yr, 
where B, in units of 10 12 G, is the surface dipole mag- 
netic field strength of a pulsar. These time scales are 
much longer than those of 10-100 yr in our model. The 
proposed test of our model with Figure [6] can be used to 
distinguish between their model and ours. 
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